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We study the statistical mechanics of quasianti- Hermitian quaternionic sys- 
tems with constant number of particles in equilibrium. We show that the explicit 
knowledge of the metric operator is not necessary for study the thermodynamic 
properties of the system. We introduce a toy model where the physically relevant 
quantities are derived. We derive the energy fluctuation from which we observe 
that for large N the relative r.m.s fluctuation in the values of E is quite neg- 
ligible. We also study the negative temperature for such systems. Finally two 
physical examples are discussed. 



1 Introduction 

Non-Hermitian Hamiltonians are currently an active field of research, motivated 
by the necessity to understand the mathematical properties of their subclasses 
, namely the pseudo-Hermitian and pseudoanti-Hermitian Hamiltonian. Also, 
to investigate the existence of a suitable similarity transformation that maps 
such Hamiltonians to an equivalent Hermitian form is important from a physi- 
cal point of view. A consistent theory of quantum mechanics demands a certain 
inner product that ensures the associated norm to be positive definite. In this 
direction there have been efforts to look for non-Hermitian Hamiltonians which 
have a real spectrum such that the accompanying dynamics is unitary. 
The interest in non-Hermitian Hamiltonians was stepped up by a conjecture of 
Bender and Boettcher [1] that PT-symmetric Hamiltonians could possess real 
bound-state eigenvalues. Subsequently, Mostafazadeh [2] showed that the con- 
cept of PT symmetry has its roots in the theory of pseudo-Hermitian operators. 
He pointed out that the reality of the spectrum is ensured [3] if the Hamiltonian 
H is Hermitian with respect to a positive-definite inner product <.,.>+ on the 
Hilbert space H in which H is acting [4] . 

In the pseudo-Hermitian representation of quantum mechanics, a quantum sys- 
tem is determined by a triplet (H,i/, 77) where H is an auxiliary Hilbert space. 
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iJ : H ^ H is a linear (Hamiltonian) operator with a real spectrum and a com- 
plete set of eigenvectors, and rj : H ^ H is a linear, positive- definite, invertible 
(metric) operator fulfilling the pseudo-Hermiticity condition.: 

H^ = ii+Hn-\ v+ = vi>o (1) 

The condition that the metric must be positive definite(quasi Hermiticity) is 
necessary for compatibility of models with the postulates of quantum mechan- 
ics [5]. 

The physical Hilbert space H phys of the system is defined as the complete ex- 
tension of the span of the eigenvectors of H endowed with the inner product 

<.,.>+:=<. I ?7. > (2) 

where < ., . > is the defining inner product of H, and the observables are 
identified with the self-adjoint operators acting in H Phys, alternatively r/-pseudo- 
Hermitian operators acting in H, [3,6]. 

The foundations of quaternionic quantum mechanics(QQM) were laid by Finkel- 
stein et.al., in the 1960's[7]. A systematic study of QQM is given in [8], which 
also contains an interesting list of open problems. 

It is worth mentioning that while in both QQM and CQM(complex quantum 
mechanics) theories, observables are associated with self-adjoint (or Hermitian) 
operators, the Hamiltonians are Hermitian in CQM, but they are anti-Hermitian 
in QQM, and the same happens for the symmetry generators, like the angular 
momentum operators. 

On the other hand as we mentioned earlier theoretical framework of CQM has 
been extended and generalized by introducing the pseudo-Hermitian operators. 
By the same motivation if one wishes to extend and generalize the theoretical 
framework of standard quaternionic Hamiltonians and symmetry generators, 
one needs to introduce and study the pseudoanti-Hermitian quaternionic oper- 
ators. 

Moreover, the theory of open quantum systems can be obtained, in many rele- 
vant physical situations, as the complex projection of quaternionic closed quan- 
tum systems [9]. 

Experimental tests on QQM were proposed by Peres [10] and carried out by 
Kaiser et.al., [11] searching for quaternionic effects manifested through non- 
commuting scattering phases when a particle crosses a pair of potential barriers. 
See also [12]. A review of the experimental status of QQM can be found in [13]. 
By definition a quaternionic linear operator H is said to be ( 77-)pseudoanti- Her- 
mitian if a linear invertible Hermitian operator r] exists such that r}Hri~^ = —H^. 
If 7] is positive definite, H is said quasianti-Hermitian. 
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2 pseudoanti-Hermitian quaternionic systems. 



A quaternion is usually expressed as q = qo + iqi + jq2 + kq^ where go, 1,2, 3 € -R 
and = j'^ = fc^ = ijk = = —ji = k, and with an involutive anti- 

automorphism(conjugation) such that, q ^ q = qo — iqi — jq2 — kqz- 
The density matrix associated with a pure state ?/> belonging to a quaternionic 
n— dimensional right Hilbert space Q" is defined by : = and is the 

same for all normalized ray representatives [14]. 

Denoting by the adjoint of an operator Q with respect to the pseudo-inner 
product (., .)t7 = (•)»?•)> (where (.,.) represent the standard quaternionic inner 
product in the space Q"). We have 

Qt = r?-iQt^ (3) 
so that for any ry-pseudo-Hermitian operator i.e., satisfying the relation 

r?Q77-i = Qt (4) 

one has, Q = Q*- These operators constitute the physical observables of the 
system. If Q is r?-pseudo-Hermitian, Eq.(4) immediately implies that r]Q is 

Hcrmitian. so that the expectation value of Q in the state | ip > with respect to 
the pseudo-inner product (., .)^ = (.,77.) can be obtained [14]: 

< ^ I r?g I V >= ReTr{\ ^ >< ^ \ rjQ) = ReTr{pQ) (5) 

where p =| V >< V' I V- 

More generally, if p denotes a generic quaternionic (Hermitian, positive defi- 
nite) density matrix, we can associate it with a generalized density matrix p by 
means of a one-to-one mapping in the following way: 

P = PV (6) 
and obtain < Q >n= ReTr{pQ)[l4]. Note that p is ?7-pseudo-Hermitian: 

P = VP = VPV~^ (7) 

Let us then consider the space of quaternionic quasi-Hermitian density matrices, 
that is the subclass of ry— pseudo-Hermitian density matrices with a positive rj. 
Thus, an (Hermitian) operator 6 exists such that rj = G^, and the r?— pseudo- 
Hermitian density matrices are positive definite. Then the inner product (., .),; = 
(., r].) introduced in the Hilbert space is positive and the usual requirements for 
a proper quantum measurement theory can be maintained. 
The most general 2-dimensional complex positive r] operator is given by : 
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where x,y G R, z £ C and xy z ^[14]. The choise of complex metric oper- 
ators can be justified as follows : it is shown in [14] that the complex projections 
of timc-dcpcndcnt 77-quasianti-Hermitian quatcrnionic Hamiltonian dynamics 
are complex stochastic dynamics in the space of complex quasi-Hermitian den- 
sity matrices if and only if a quasistationarity condition is fulfilled, i.e., if and 
only if rj is an Hermitian positive time- independent complex operator. 
Now let H describes an ensemble of a huge but fixed number of independent 
particles. In what follows we study the statistical description of the system in 
equilibrium. The state of the system is characterized by a density matrix p 
which can be written in energy representation as 

P = ^W„|n)(n|, Trp < 00 (8) 

n 

In equilibrium, the density matrix operator is solution of Bloch equation: 

with initial condition p(0) = 1 and /? = -^jy , k is the familiar Boltzmann constant. 
Its formal solution p = e~^^ fixes the coefficients to be W„ = exp{—(3En). 
Normalization factor of the density matrix depends on the inverse temperature 
/? and is called partition function. 

Z = Trp (10) 

It plays an important role in thermodynamics of the system because it allows 
direct computation of thermodynamic quantities. 

Using the definition of the expectation value of an operator in quasianti- Hermitian 
picture i.e., Eq.(5) one can show that : 

Z = J2{n\p\n) = ^{n\p\n)=Z (11) 



where \n), pand Z are energy eigenkets, density operator and partition func- 
tion in quasianti-Hermitian quaternionic picture respectively. This means that 
the partition function is the same in both cases. So we can derive the thermo- 
dynamic properties of a physical system in equilibrium in pseudoanti-Hermitian 
quaternionic picture without explicit knowledge of the metric operator O. 
On the other hand Bloch equation (9) resembles Schrodinger equation for evo- 
lution operator in QQM. The link can be established by the substitution (3 ^ t. 
Working in quasianti-Hermitian picture the equation (9) reads: 

^ = -HU. U{0) = 1 (12) 
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where we defined U(t) = p(t). In quasi-Hermitian picture. In x-representation, 
we have 

p{xi,X2) \f3^t= U{xi,X2\t) =< xi I e~"* I X2 >=: G{xi,X2;t). (13) 

the propagator (13) is of fundamental importance because it allows to com- 
pute partition function [Z = J G{xi,X2;—p)dx) with use of standard tech- 
niques. The analogous equations in the case of pseudo-Hermitian Hamiltonian 
have been derived in [15]. 



3 Statistical mechanics of quasianti-Hermitian 
quaternionic systems. A Toy model. 

We consider an ensemble of systems each consisting of N distinguishable par- 
ticles without mutual interaction. The partition function can be written as 
follows : 

Z=iZif (14) 

where Zi may be regarded as the partition function of a single particle in the 
system. In the quasianti-Hermitian quaternionic picture, the subsystem is de- 
scribed by the following most general Hamiltonian : 




Here a,b and c are three arbitrary quaternion. The requirement that the 
Hamiltonian be pseudoanti-Hermitian quaternionic with respect to rj i.e. riHri~^ = 
—H^ gives : a* = —a, b* = —b and d = —^c*, where t] is given by : 

where we have chosen z = 0, x^ = a and t/^ = 7 in the matrix rj introduced 
in section 2. The energy levels of the system may depend on the external 
parameters e.g. on the volume in which the particle is confined. 
We shall find the partition function by solving the Eq.(12). By separating the 
Hamiltonian into its diagonal and non-diagonal part we have 

H = Ho + H' (15) 

where : 

( b)' -L* ) 
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Let us introduce the following transformation : 

Uo{t) = exp{-Hot) (16) 
we also define a new operator in the interaction picture as follows: 

Ui{t) = u^^mmoit) (17) 

Then Eq.(12) changes into the following form : 

= -H'j{t)Ur{t), urn = 1 (18) 

where Hj{t) = Uo{ty H'Uo{t). One can find the solution iteratively by integrat- 
ing both sides of the equation. 

Ui{t) = 1- f dt'Hi{t') + /* dt' f dt" H'j{t')H'j{t") + ... (19) 
Jo Jo Jo 

where we have neglected the higher order terms, then we have : 



/ l + --S;i--r^(e^''"'''*-l) « / (a-6)t _ 

' ^ \ d (e-(«-f)* - 1) 1 _ _ « cl^/g-(a-6)t _ 1) 

\ a a— ^ ^ 7 a— 7 (a— o)^ ^ ' 



The approximate value of the partition function is as follows : 

Zi(/3) = TT{U^{t)Ui{t)) |t^;3~ e-'^^ + e-"^ - '^^Me-^P - e-l') (20) 

7 a — 

Now the calculation of the thermodynamic quantities is straightforward. For 
the Helmholtz free energy of the system we have : 



A = -iVfcT In Zi = - ^ InZi = Infe"'^'^ + e"'"' - -^^ifT^^ - 6""'^)] 
p p 'y a—b 

(21) 

The entropy of the system is : 

S = 7Vfcln[e-"^ + e-''^ - ^-£^l3{e-^^ - e""'^)] 

\a{a -h)- ^c^ + ^c^Pa]eJ^'' + \b{a - 6) + ^c^ - ^c^me'^'' 

[(a-6) + ^c2/3]e/3'' + [(a-6)-^c2/3]e/3'» ^ ' 



6 



(24) 



The internal energy of the system is given by : 

e^^[a{a - h + f c^^) - f c^] + e'3'^[5(a - 6 - f c^/?) + f c^] 

Tj = /V 2 2 2 2 

e^b^a - 6 + ^c2/3] + e/3'^[a - b - ^c^(3] ^ ' 

The specific heat per particle Cy which describes how the temperature 
changes when the heat is absorbed while volume V of the system remains un- 
changed is given by : 

n — (dU \ _ -0^k (dU \ _ 

p2j^ _e/3(«+b){(a _ 6)2 [(a _ 5)2 _ (^c2/3)2 - 4^c^] + 2{^c^f} + {'^^(?f\e^&<>- + e^<^^\ 
~N {el^b^a -b+ ^c^j3] + e/5«[a - 6 - ^c2/3]}2 

The pressure of the system is as follows : 

p__/OA^ _ AT + f c'b'+2f c(a-b)c'] 



\-{a - bfh' - ^c^il - (iia - b))b' + ^c^a' - 2^c{a - 6)c'] , 
[{a - 6)2 + ^c2/3(o _ 6)]e-a/3 + [(a - 6)2 - ^c2/?(a - 6)]e-''/3 ^ ^ 

where the prime symbol represents the derivative with respect to the volume, 
" dv 

4 Energy fluctuation. 

We first write down the expression for the mean energy : 
U=<E>= ^rEr9reM-PE) 

J2r9rexp{-(3Er) 

where is the multiplicity of a particular energy level E,.. By differentiating of 
the expression of the mean energy with respect to the parameter f3, we obtain : 

^ = -{{E^) - {Ef} (27) 

whence it follows that : 

{{AEf) = (E^) - {Ef = = fcT^(^) = kT^Cy (28) 
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So, for the relative root-mean-square fluetuation in E, we have : 

Vii^m ^ VkT^ ^ 1 (3^kf{a,b,c,a,P,j) ^ J_ 
(E) U y/N g{a,b,c,a,f3,j) y/N 

where f{a,b,c,a,P,j) is the numerator of Eq.(24) and g{a,b,c,a, is the 
numerator of Eq.(23). As we observe it is 0(-/V~2), N being the number of 
particles in the system. Consequently, for large 7V(which is true for every sta- 
tistical system), the relative r.m.s fluctuation in the values of E is quite neg- 
ligible. Thus for all practical purposes, a system in the canonical ensemble in 
pseudoanti-Hermitian quaternionic picture, has an energy equal to or almost 
equal to the mean energy U ; the situation is therefore practically the same as 
in the micro canonical ensemble. 



5 Negative temperature. 

Let us consider our system from the combinatorial point of view. The question 
then arises: in how many different ways, can our system attain a state of energy 
E?. This can be tackled in precisely the same way as the problem of the random 
walk. Let be the number of particles with energy E+ and N- with energy 
E-; then 

E = E+N++E-N-, N = N+ + N- (30) 

Solving for A^_|_ and , wc obtain : 

E-NE. E-NE+ , , 

Ar+ = -= N. = -± (31) 

The desired number of ways is then given by the expression: 

n{N, E) = = , E-NE+ .u E-NE_ ., (32) 

whence we obtain for the entropy of the system : 

S{N,E) = k\nn ~ 

,r.., .r E-NE+^ ,E-NE+, E-NE^^ ,E-NE^,^ 

^[^^^^ - e7-e:'''^^7^^ - 1^(^73^)] 

The temperature of the system is then given by : 

1 .5-5, k , , E-NE_^ 

f = (a^^- = et^'^'^-e^ne;^ 

from Eq.34,we note that so long as ii^ < ^(^-i- +^-)) T > . However the same 
equation tells us that ii E > ^{E^ + E_), then T < 0. Lets examine the matter 
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a little more closely. For this purpose we consider as well the variation of the 
entropy S with the energy E. We note that for E = NE- , both S and T vanish. 
As E increase, they too increase until wc reach the special situation where 
E = +-£'-)• The entropy is then seen to have attained its maximum value 

5 = Nkln2, while the temperature has reached an infinite value. Throughout 
this range, the entropy was a monotonically increasing function of energy,so T 
was positive. Now as E equals [y{E+ + E_)]^, (^) becomes 0_ and T becomes 
— oo. With a further increase in the value of E, the entropy monotonically 
decreases; as a result, the temperature continues to be negative, though its 
magnitude steading decreases. Finally, we reach the largest value of E, namely 
A''£+,where the entropy is once again zero and T = — 0. 

6 Physical examples. 

6.1 A spin one half system in a constant quasianti-Hermitian 
quaternionic potential. 

We now consider a two-level quantum system with a quasianti-Hermitian quater- 
nionic Hamiltonian H — Ha+jHr;. denotes the free complex anti-Hcrmitian 
Hamiltonian describing a spin half particle in a constant magnetic field[14]. 




and jHfj is a purely quasianti-Hermitian quaternionic constant potential : 



■' \ jvx J 

We note that H = Ha -|- is ry'-quasianti-Hermitian quaternionic i.e., 
Tj' Hri'~ = —H\ where : 




The eigenvalues and the corresponding biorthonormal eigenbasis of the quater- 
nionic Hamiltonian H are[16]: 
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iE± = ± v) (35) 



and : 



We note that the metric r]' is a special case of our metric and can be obtained 
by substituting a = and 7=1. We also note that the Hamiltonian H is 
a special case of our general Hamiltonian and can be obtained by substituting 
a = i^,b = —i^,c = and d = —^c* . So the thermodynamic properties 
of this system can be obtained straightforwardly from equations (21-25), for 
instance we have : 

S = A^fcln[2(cosh ^ - ^ sinh ^)] 

sij^ + 2v^ sinh + (3ujv'^ cosh , ^ 

+Nk(} 2 — 2 — t: — 2_ 3g 

2a; cosh -2/3t;2 sinh ^ ^ ' 



^uP' sin + 'Iv- riiuli + s'lxv- cosh 4r 



U = iV 2 ^ 2 ^ ^ 2_ C37\ 

2a;cosh^i^-2/3t;2sinhif 



Moreover the negative temperature discussed in previous section is applicable 
to this system which we study now. 

Let Nj^ be the number of particles with energy = ^ + and N- with energy 
E_ = ^ — then we have : 

E = E+N+ + E_N_, N = N+ + N_ (38) 

The number of ways this system can attain a state of energy E is given by : 

N\ Nl 

" NJnJ- " / E-N(f+v) .,, E-N(f-v) ., (39) 

whence we obtain for the entropy of the system : 



S{N,E)=k\nn 

E-N(%+v) , E-N(%+v)^ E-N(%-v) E- -v)^, 
k[N In N+ !^ ^In — -) — ^In — -)] 

(40) 



The temperature of the system is then given by : 

1 ^^^^ ^ 1 r E-m^-v) 
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We note that so long as ^ < f w, T > . However if ^ > ^co, then T < 0. 
Now we consider as well the variation of the entropy S with the energy E. We 
note that for E = N{j — v), both S and T vanish. As E increase, they too 
increase until we reach the special situation where E = ^lu. The entropy is 
then seen to have attained its maximum value S = Nk In 2, while the temper- 
ature has reached an infinite value. Throughout this range, the entropy was a 
monotonically increasing function of energy,so T was positive. Now as E equals 
becomes 0_ and T becomes — oo. With a further increase in the 
value of E, the entropy monotonically decreases; as a result, the temperature 
continues to be negative, though its magnitude steading decreases. Finally, we 
reach the largest value of E, namely + w),where the entropy is once again 
zero and T = —0. 

As we mentioned earlier the physical system is a spin half particle in the 

quatcrnionic potential. How to produce in laboratory quatcrnionic potentials 
has been investigated by Peres, Mc Kellar, Brumby and others [10,11,12,17,18]. 
This system is the quasianti-Hermitian quaternioinc analogue of the ordinary (Hermitian) 
magnetic system i.e., a system composed of spin half particles in a constant mag- 
netic field which is a famous system showing negative temperature discussed in 
the text books, see e.g. ,[20]. Each dipole has a choise of two energy = —fiH 
and £J+ = +fiH. 



6.2 Optimal entanglement generation. 

The next system we will study is composed of two C qubits A and B, parametrized 
by the Bloch vectors a = (ai, 02, 0,3) and b = (61, 62, 63), respectively. The most 
general Hamiltonian for two qubits (in the interaction picture) can be written 
as: 

H = Y.Ci<^f(S}<^f (42) 

i 

where the parameters Ci , C2 and (3 are constant if there is no free evolution for 
indivitual qubits. Assume now we have initially two complex pure states, with 
Bloch vectors a = (1, 0, 0) and b = (0, 1, 0), respectively. We also assume some 
interaction associated with the following parameters : 

C3 = l,Ci=C2 = (43) 

It is shown in [19] that the corresponding anti-Hermitian quatcrnionic Hamil- 
tonian has the following form: 
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It is easy to check that the Hamiltonian is anti-Hermitian quaternionic with 
respect r] given in section 3. 

Again we note that this Hamihonian is a special case of our general Hamiltonian 
and can be obtained by substituting a = — 2je~*'^, 6 = and c = 0, so its ther- 
modynamic properties can be obtained from equations (21-25) by substitution 
the values of a, h and c. 

We know that any two- level system can be used as a qubit. Several physical 
implementations which approximate two-level systems to various degrees were 
successfully realized. The following are three examples of physical implementa- 
tions of qubits(the choices of basis are by convention only). 

1) . Optical lattices (Atomic spin), |0): up ; |1): down. 

2) . Joscphson junction(Superconducting charge qubit), |0): Uncharged super- 
conducting island (Q=0) ; |1): Charged superconducting island (Q=2e, one 
extra Cooper pair). 

3) . Quantum dot(Dot spin), |0): up ; |1): down. 



7 Conclusion. 

We have studied the statistical mechanics of quasianti-Hermitian quaternionic 
systems in equilibrium in the framework of canonical ensemble. We derive all 
the physically relevant quantities without explicit calculation of the metric op- 
erator or the spectrum of the Hamiltonian. We derive the energy fluctuation 
from which we observe that for large N, the relative r.m.s fluctuation in the 
values of E is quite negligible. We also study the negative temperature for such 
systems. Two physical examples are discussed. 



8 Acknowledgment. 

I would like to thank A. Tatar for his kind assistance. 



9 References. 

1. CM. Bender, S. Boettcher, Phys. Rev. Lett.80(1998)4243. 

2. A. Mostafazadch, J. Math. Phys.43(2002)205. 

3. A. Mostafazadch, A. Batal, J. Phys. A37(2004)11645, J. Phys. A38(2005)3213. 

4. B. Bagchi et.al., J. Phys. A38(2005)L647. 

5. M. Znojil, H. B. Gcycr, Phys. Lett. B640(2006)52. 

6. A. Mostafazadeh, Phys.Rev. D 76, 067701 (2007). 

7. D. Finkelstein et.al., J. Math. Phys.3(1962)207; 28(1963)788. 

8. S. L. Adler, Quaternionic quantum mechanics and quantum fields, Oxford, 



12 



New York, 1995, and references therein. 

9. M. Asorey and G. Scolarici, J. Phys. A 39 (2006) 9727; M. Asorey, G. Sco- 
larici and L. Solombrino, Theor. Math. Phys. 151 (2007) 733; M. Asorey, G. 
Scolarici and L. Solombrino, Phys. Rev. A 76 (2007) 12111; A. Kossakowski. 
Rep. Math. Phys. 46 (2000) 393; A. Blasi et.al., quant-ph/0407158. 

10. A. Peres, Phys. Rev. Lett. 42(1979)683. 

11. H. Kaiser et.al., Phys. Rev. A29 (1984)2276. 

12. S. L. Adler, Phys. Rev. D37 (1988)3654. 

13. S. P. Brumby, G. C. Joshi, Chaos, Solitons Fractals 7(1996)747. 

14. G. Scolarici, SIGMA3(2007)088, G. Scolarici, L. Solombrino, arXiv:0711.1244. 

15. V. Jakubsky, Mod. Phys. Lett. A22(2007)1075. 

16. G. Scolarici, J. Phys. A:Math.Gen.35(2002)7493. 

17. A. J. Davies and B. H. J. McKellar, Phys. Rev. A40(1989)4209; 46(1992)3671; 
A. J. Davies, Phys. Rev. D41 (1990)2628. 

18. S. P. Brumby, et.al., Phys. Rev. A51(1995)976. 

19. Asorey, et.al., Phys. Rev. A76(2007)012111. 

20. R. K. Pathria, Statistical Mechanics, Pergamon press, 1998. 



13 



